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Abstract. In this paper techniques of twist map theory are applied to the annulus 
maps arising from dual billiards on a strictly convex closed curve T in the plane. It 
is shown that there do not exist invariant circles near T when there is a point on F 
where the radius of curvature vanishes or is discontinuous. In addition, when the 
radius of curvature is not C 1 there are examples with orbits that converge to a point 
of r. If the derivative of the radius of curvature is bounded, such orbits cannot exist. 
The final section of the paper concerns an impact oscillator whose dynamics are the 
same as a dual billiards map. The appendix is a remark on the connection of the 
inverse problems for invariant circles in billiards and dual billiards. 

§0 Introduction 

Dual billiards is a dynamical system defined on the exterior of an oriented convex closed 
curve T in the plane. If z is a point in the unbounded component of M 2 — T, then its image 
under the dual billiards map $ is the reflection about the point of tangency in the oriented 
supporting line to T (see Figure 0.1). It is clear that <3> is an area-preserving, and if T is 
strictly convex, it is homeomorphism. Further, $ is affine invariant, i.e. transforming the 
convex curve by an affine map simply transforms the entire homeomorphism. If the initial 
curve is an ellipse, then its exterior is foliated by concentric homothetic ellipses that are 
$- invariant. 

c5(z) 




Figure 0.1: The dual billiards map. 

The invention of dual billiards is credited to B. H. Newman in [Ms2]. The nomencla- 
ture outer billiards or cobilliards is also used. The name dual billiards is apt because this 
dynamical system is in certain ways the dual of the usual billiards. Indeed, many of the 
theorems given here are analogs of theorems for billiards. There are, however, fundamental 
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differences between the systems. Perhaps the most important is that the phase space for 
the dual billiards map is non-compact and has infinite area. An obvious obstruction to any 
explicit duality between the two systems is that the billiards map defined using an ellipse 
has period-two points of elliptic and saddle type. 

One of the goals of this paper is to show how dual billiards fits into recent developments 
in the theory of area-preserving twist maps of the annulus due mainly to Mather. Once it 
is known that a dynamical system falls into this class there are a wealth of results available 
which yield a great deal of information about the dynamics of the map. There are technical 
differences, however, between the twist maps that arise in dual billiards and the usual classes 
of twist maps that are studied in the literature. 

Section 1 begins by specifying a class of twist maps called half-cylinder twist maps that 
includes those from dual billiards. The rest of §1 is devoted to describing some standard 
tools and results of twist map theory and remarking on how these can be adapted to the 
new class. We only describe that part of twist map theory that is relevant. For a broader 
perspective the reader is referred to the various papers of Mather (summarized in [M5]), as 
well as [Bg], [Ms3], [Ms4], chapter 1 of [H] and chapter 10 of [MH]. There are surveys with 
a more physical point of view in [Mc] and [Me]. 

Section 2 give necessary background material on convex curves and in §3 it is confirmed 
that, in the coordinates introduced in §2, the dual billiards map is indeed a half-cylinder 
twist map. The next section deals with topological circles that are invariant under the dual 
billiards map. The question of the existence of invariant circles is central in twist map 
theory. An invariant circle provides a barrier through which orbits cannot pass, and thus 
the existence of invariant circles is intimately connected with stability questions. Section 
4 begins with two results of R. Douady that give existence of invariant circles near T and 
near infinity. These results require a certain smoothness as they ultimately depend on KAM 
theory. 

The inverse problem for invariant circles is considered next: given a convex curve Ti is 
there a T whose dual billiards map has Ti as an invariant circle. The solution to this problem 
is given by a geometric construction called the area envelope. The analog of this construction 
for billiards is discussed in the appendix where it is pointed out that both the billiard and 
dual billiard constructions are a consequence of the constancy of an area function (Mather's 
big "H") that can be associated with an invariant circle of a twist map (cf. Proposition 1.3). 

The next results in §4 concern the non-existence of invariant circles near T when there 
is a point on T where the radius of curvature vanishes or is discontinuous. In addition, when 
the radius of curvature is not C 1 there may also exist orbits that converge to a point of T. 
These so-called crash solutions are shown to exist for certain examples in §5. It is also shown 
that if the derivative of the radius of curvature is bounded, no such orbits exist. 

There are no known examples of strictly convex curves for which the dual billiards map 
lacks invariant circles near infinity. This problem (raised by Moser) remains one of the most 
important outstanding problems in the theory of dual billiards. 

The final section of the paper concerns an impact oscillator whose dynamics are the 
same as a dual billiards map. Consider a wall that is periodically moving with position 
given by p(t) where pit) satisfies p + p = p(t) for some 27r-periodic forcing function p with 
/ exp(it)p(t) = 0. A simple harmonic oscillator with position specified by x(t) satisfying 
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p(t) 



x(t) 

Figure 0.2: The impact oscillator. 

x + x = is assumed to have perfectly elastic collisions with the wall. Figure 0.2 shows one 
possible (idealized) realization of the system. On the left is a fixed wall whose position is 
set at the origin. A pair of pistons periodically drive a second wall that is to the right of 
the fixed wall. A flat plate is connected to the fixed wall by a spring that passes through 
a hole in the moving wall. We monitor the motions of the plate as it moves away, collides 
with the moving wall, moves away, etc. The return map to collisions with the moving wall 
in appropriate coordinates turns out to be identical to dual billiards on a convex curve T 
whose radius of curvature function p is the same as the forcing function for the wall. 

As a consequence of the equivalence of the two dynamical systems one can use results 
about dual billiards to provide information about the oscillator. For example, as a conse- 
quence of the Birkhoff-Mather stability theorem (see Remark 4.3.2) when there do not exist 
^-invariant circles near the curve T, there must exist orbits that converge to T in forward 
and backward time. Thus as a consequence of Theorem 4.3, if there is an instant at which 
the forcing function is zero or discontinuous, then there are solutions for the impact oscillator 
that converge to the wall. In addition, it follows from Theorem 5.1 that there are impact 
oscillators that have crash solutions that converge to the wall in finite time. These solutions 
do not exist if the forcing function is C 1 (Theorem 5.2). 

The situation near infinity is perhaps of more physical interest. As a consequence of 
Theorem 4.1, for sufficiently smooth convex curves the dual billiards map has invariant 
circles near infinity. Thus for sufficiently smooth forcing functions the impact oscillator has 
no unbounded solutions. Because of the lack of unbounded solutions these oscillators are 
called stable. There have been numerous other impact oscillators studied in the literature. 
Their stability depends on the details of the oscillator, eg. on the restoring force of the 
particle to the wall, the relation of the frequency of the forcing to that of the free oscillator, 
etc. Again using the Birkhoff-Mather stability theorem, the non-existence of invariant circles 
near infinity for dual billiards on some convex curve would imply unbounded solutions for 
the corresponding impact oscillator. This gives a physical motivation for the question of 
Moser noted above. 

This paper has several purposes, and it contains material of different types. A first goal 
is to survey known results from the point of view adapted here. For this reason there is an 
overlap in basic material with other papers that have appeared on dual billiards, for example, 
[GK], [Tl], [T2], and [T3]. A second goal is to give a mathematical framework for future 
work on dual billiards. Thus there are somewhat technical sections (eg. §1.1) that could be 
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skipped by the reader with more physical interests. Despite these rigorous goals we also want 
to preserve the essentially simple and intuitive geometric character of the problem. Thus 
there also is descriptive material. Much of this material also serves to illustrate geometric 
heuristics associated with twist map theory. 

Acknowledgments: Most of the results in this paper were obtained while the author was a 
Postdoctoral Research Fellow at the Mathematical Sciences Research Institute in Berkeley, 
1988-89. He would like to thank them (belatedly) for there support. Thanks also to Danny 
Goroff, Dick Hall, and Robert MacKay for useful and stimulating conversations. 

§1 Half-cylinder twist maps 

In this section we introduce a class of area-preserving monotone twist maps called half- 
cylinder twist maps, and point out that the standard results from twist map theory apply 
to this class. We will see in §3 that a dual billiards map in the appropriate coordinates is a 
half-cylinder twist map. 

§1.1 Definition of half-cylinder monotone twist maps. Let the one point compactifi- 
cation of the ray [0, oo) be denoted [0, oo] with the added point being labeled oo. Let C be 
the open half-cylinder S 1 x (0, oo) and C* be S 1 x [0, oo], with the circle S 1 = M/(2ttZ). The 
space C* is homeomorphic to the compact annulus, but we maintain the infinite measure. 
The universal covers of these spaces are denoted by C and C* , respectively. A tilde will 
always indicate the lift of a point, function, set, etc. to the universal cover. The map ix\ 
in various contexts means the projection of a product space onto its first component. We 
will use the coordinates (#,7) for C and C* and (x, 7) for C and C*. If g is a real valued 
function, gi denotes the derivative with respect to the i th variable. 

The definition of a half-cylinder monotone twist map requires two pieces of data. The 
first is a pair of functions bo, boo '■ K — > R called the boundary maps. It is required that 60 
be nondecreasing and continuous from the right, boo be increasing and continuous, bo < &oo> 
and bi(x + 1) = bi(x) + 1 for i — 0, 00. Let V = {(x,x') E R 2 : b (x) < x' < boo{x)}. For 
i = 0, 00, A, C IR 2 denotes the graph of b{ as a function of x. Note that Ao U Aqo is contained 
in the frontier of V. 

The second piece of data needed to define a half-cylinder twist map is a C 1 -function 
h : V — > R called a generating function that satisfies: 

(Gl) (periodicity) h(x + 2vr, x' + 2vr) = h(x, x') 

(G2) (conditions on the derivatives) h\ < and fi2 > 0, and for each e > 0, the derivative 
Dh is Lipschitz on each T) e = {(x,x') E M 2 : bo(x) < x' < boo(x) — e}. Further, the mixed 
partial derivative h\2 = ^21 exists, is continuous and satisfies h\2 < 0. 

(G3) (limit behavior near Aq) For i = 1,2, 



lim hi(x, x) = 

(x,x')->Aq 



lim h(x,x'). 
(x,x')->Ao 



(G4) (limit behavior near A^) The first derivatives must satisfy 



lim li2(x,x') 



= 00 = 



lim h\(x,x') 

x->b^J-(x) 
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and for each x, 



h,2(u,v) + hi(v,w) — > as (u,v,w) — > (^(rr), rr, &oo(x))- 



Definition: A map f : C* 

f : C* — > C* i/ia£ satisfies 



C* is called a half-cylinder monotone twist map if it has a 

(1.1) 



7 = -/ii(z,a/) 
y = h 2 {x,x') 



where (x', 7') = /(x, 7) /or (x, 7) G C ; and /i satisfies properties (Gl) - ('C^j afrove 6j 
egna/ to / restricted to Rx {i}. 

As a consequence of the properties (Gl) - (G4), when restricted to C, a half-cylinder 
monotone twist map / will always be an area-preserving homeomorphism, and both / and 
f^ 1 will be locally Lipschitz. Conversely, given bi and h satisfying (Gl) - (G4), (1.1) 
implicitly defines a half-cylinder monotone twist map /. Using (1.1), Df(x, 7) is given by 



J_ ( -/ill "I 

hl2 \ h 12 ~ ^11^22 ~/l22 



1.2) 



where each of the derivatives hij is evaluated at (x,x') with x' = iri(f(x, 7)). 
Remarks 

1.1.1 The name "twist map" is given because the image of a vertical arc is always the 
graph of a function defined on a subset of the circle. From (1.2) one sees that this is a 
consequence of < 0. The generating function h can be given a geometric interpretation 
as in Figure 1.1 ([K]). The area of the region bounded by the x axis, the vertical arc above 
x', and the image under / of the vertical arc above x is equal to h(x,x'). 




h( x, x') 



Figure 1.1: The geometric meaning of the generating 
function. 

1.1.2 There are a few differences between half-cylinder twist maps and the classes of twist 
maps that occur in the literature. These differences are mainly "technical" in the sense that 
standard arguments work for the class with straightforward alterations. The consideration of 
these technical differences are forced on us by the twist maps that come from dual billiards. 
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A twist map is usually denned either on the compact annulus or on the infinite cylinder. 
In the later case one requires infinite wrapping, i. e. the image of a vertical line wraps infinitely 
around the cylinder. In contrast, half-cylinder twist map are defined on the half-infinite 
cylinder. They extend to a perhaps non-continuous map on the lower boundary, and can be 
extended to a homeomorphism on an upper circular boundary that is added at infinity. Thus 
they do not have infinite wrapping; they can be thought of as being defined on a compact 
annulus with infinite 

If a twist map / is differentiable at a point zq the positive twist at zq, denoted twist + (zq), 
is the angle (measured clockwise) from a vertical vector based at f(zo) to the image under 
Df of a vertical vector based at zq. The lift of a half-cylinder twist map / may not be 
differentiable, but the image of the vertical line through zq = (a?o,7o) * s always the graph 
of the Lipschitz function of x', h,2(xQ,x'). In this case we let twist + (zo) = aicot(K) where 
K is the best Lipschitz constant for this function at x'q = ni(f(zo)). The negative twist 
at zq, denoted twist" (zq), is the twist at zq of S o f where S(x, r y) = (— x, 7). The twist 
at zq is twist(zo) = min{twist + (zo),twist-(zQ)}. Note that if h is twice differentiable, 
twist(zo) = min{arcot(/in(xo, xq), arcot(/i22(^o> x'q)}, where xq = 7ri(/ _1 (,2o)). 

In twist map theory it is usually assumed that the twist is bounded away from zero. For 
half-cylinder twist maps the condition given in (G2) requiring that Dh is Lipschitz in T> e 
insures that the twist is bounded away from zero in neighborhoods of the lower boundary of 
C. However, / itself may only be only locally Lipschitz as its Lipschitz constant can go to 
infinity as one approaches the lower boundary. In a half-cylinder twist map the twist will 
not be bounded away from zero as one approaches infinity as the map extends to the circle 
at infinity. 

§1.2 The variational formulation and dynamics of twist maps. Note that (1.1) 
says that h is a generating function of a canonical transformation in the sense of classical 
mechanics (eg. [MH], page 47). The function h plays another role and is sometimes called the 
action as it can be used in a discrete variational formulation. What follows is the rudiments 
of this theory adapted for half-cylinder twist maps. For a full treatment the reader is referred 
to [Bg]. 

A configuration is an element x = (xi)i £ z of the set M z of bi-infinite sequences of real 
numbers (the nomenclature "configuration" comes from solid state physics). Given a gener- 
ating function h and boundary maps b% as above, let X = {x £ R : (xi, Xi+\ ) £ £> for all i}. 
Elements of X are called allowable configurations. The generating function h can be extended 
to finite segments of allowable configurations via 

fe-i 

h(xj, ... ,x k ) = J2 h ( x h Xi+i). 
i=j 

An allowable configuration x is called stationary if for all pairs (j, k) with j < k, h(xj, . . . , x^) 
is stationary with respect to variations with fixed endpoints. The configuration is called 
minimizing if for all pairs (j, k) with j < k, h(xj, . . . , x k ) is minimal with respect to variations 
with fixed endpoints. Minimizing configurations are always stationary. A configuration is 
stationary if and only if 

h,2(xi-i,Xi) + h\(xi,Xi + \) = (1.3) 
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for all i, or as a consequence of (1.1), precisely when a; is a sequence of x coordinates of an 
orbit of the map /, i.e. for all i, X{ = ni(f(xo, 70)) for some initial (xo,7o) G C. 

A configuration is monotone if for each fixed (m, n) G 1? either Xi+ m > rcj + n for all i, or 
else, Xj+ m < Xi + n for all i. A basic featureof the theory is that minimizing configurations 
are always monotone. An orbit (xi,7i) of / is monotone if the sequence (xj) is. An orbit 
(0i, 7i) of / is monotone if it has a lift to C that is monotone. Roughly speaking, a monotone 
orbit is one on which / preserves the radial order. An invariant set is called monotone if 
each of its orbits is. 

Given the lift of a half-cylinder twist map / the rotation number of a point z G C* is 

rot{i jy. = 

v 27m 

if the limit exists. The rotation number of a point z G C* is the rotation number of one 
(and thus all) of its lifts. Note that this is defined only up to a choice of lift of /, i.e. up 
to integer translation. The rotation set of / is rot(f) := {rot(z, /) : z G C*}. The rotation 
number of a point under a circle map is defined similarly, only in this case, if g : S 1 — > S 1 is 
nondecreasing (but not necessarily continuous), then the rotation number exists and is the 
same for all points in the circle. This common number is called the rotation number of the 
map, and is denoted rot(g). Another basic feature of twist map theory is that monotone 
orbits always have a rotation number. This is basically because they behave like orbits of 
circle maps. Note that any nondecreasing G : R — > R that satisfies G(x + 2tt) = G(x) + 1-n 
may also be assigned a rotation number that will be denoted rot(G). 

A Denjoy minimal set for a half-cylinder twist map is a compact, invariant set X so that 
/ restricted to X is conjugate to the exceptional minimal set in a Denjoy counterexample on 
the circle (see eg. [Dv], page 108). The phrase "invariant circle" in this paper always means 
a homotopically nontrivial circle that is invariant under the map. Since / restricted to an 
invariant circle, Q, is a circle homeomorphism, it has a rotation number that will be denoted 
rot(Q). 

Theorem 1.1: If f is a half cylinder twist map with boundary functions hi, then: 

(a) (Aubry- Mather) The rotation set of f is rot(f) = [rot (bo), rot (boo)]. For eachp/q G 
rot(f) with p and q relatively prime, there is a monotone periodic orbit with rotation number 
p/q and period q. For each u G rot(f) — Q, there is an invariant circle or Denjoy minimal 
set Xcj that is monotone and each x G X^ has rotation number uj. 

(b) (Birkhoff) IfQ C C is a f -invariant circle, then Q is the graph of a Lip schitz function 
S 1 — > (0,oo). Further, if C denotes the set of invariant circles (including S 1 x {0}) with the 
Hausdorff topology, then C is closed and the map rot : C — > R is continuous and monotonic 
in the sense that Q2 being contained in the open annulus bounded by fli and S 1 x {0} implies 
rot(Qi) > rot{Vl2) > rot(bo). 

(c) (Mather) If f has a invariant circle Vt C C, then for almost every (#0,70) e ^» 

h\i{x Q ,x\) + h 2 2(x-i,x ) > 

where (xj,7j) = (rr ,7o). 
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Proof: Part (a) can be proven using variational methods in which case the monotone 
orbits are obtained as minimizing orbits (eg. [Bg]). It may also be proven using topological 
methods (eg. [MH]), or using monotone recursion relations ([A]). In part (b), the statement 
that invariant circles are Lipschitz graphs can be proved as in [Ml] and [M2], or chapter 1 
of [H]. The statement that the set of invariant circles is closed in the Hausdorff topology 
requires the observation (explained in Remark 1.2.2 below) that condition (G2) implies a 
uniform bound on the Lipschitz constant of functions u : S 1 — > (0, oo) whose graphs are 
invariant circles in S 1 x [0, K\. The monotonicity of the rotation numbers on the circles is a 
standard consequence of the twist hypothesis. 

To prove (c), let a denote the homeomorphism of the real line that is first component 
of / restricted to the lift of the invariant circle, i.e. a(x) = iri(f(x,u(x))). Since / is locally 
Lipschitz and u is Lipschitz, a and a -1 are Lipschitz, and thus their derivatives exist almost 
everywhere. Equation (1.2) yields that 

hi(x, a(x)) + h 2 (a~ 1 (x),x) = 0. (1.4) 

Part (c) then follows by differentiating and using hn < 0. □ 
Remarks: 

1.2.1 A geometric interpretation of the condition in Theorem 1.1(c) can be given in 
terms of "tumbling tangents" (cf. [Mc] and [MP]). If h is twice differentiable, using (1.2) one 
sees that the slope of the image of the vertical tangent vector at the point (x, 7) is given 
by h22(x,x'), where x' = 7T\(f(x, 7)). Similarly, under f~ 1 , the slope of the image of the 
vertical tangent vector at the point (a/, 7') is — hn(x, x'). 

Now assume there is an invariant circle that passes through the three points (#i,7i) for 
i = —1,0, 1 where (xj+i, 7^+1) = f(xi,ji). Theorem 1.1(b) says that the circle is the graph 
of a Lipschitz function, and thus the unit tangent vector in the direction of the circle exists 
almost everywhere; let us assume that it exists at our triple. Since the circle is invariant, 
the collection of these circle directions is invariant under the induced action of Df on the 
unit tangent bundle. Now since / is a twist map, the image of a unit vertical vector always 
points in the positive x direction. Because the three points are on an invariant circle and the 
circle directions are invariant, the second iterate of the vertical vector at (x_i,7_i) cannot 
rotate past the downward vertical. 




Figure 1.2: The tumbling tangents criterion for the 
non-existence of invariant circles. 
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Rather than apply these considerations to the second iterate, it is usually more convenient 
to iterate forward from (x~i, 7-1) and backwards from (x\, 71). Since the slope of the forward 
image of a vertical vector is h\\ and a backward image is — /122, the "no tumbling tangents" 
condition is exactly that given in Theorem 1.1(c). Figure 1.2 shows an example in which 
there can be no invariant circle containing the three points. 

By using more iterates, these considerations can be refined to obtain the converse KAM 
theory of MacKay and Percival ([MP]). The macroscopic version of tumbling tangents is the 
existence of non-monotone periodic orbits as in [B+H]. 

1.2.2 The regularity conditions imposed on half-cylinder twist maps were chosen so that 
the standard proofs of Theorem 1.1(b) and (c) would go through with little change. In order 
to differentiate (1.4) almost everywhere, we need that / is locally Lipschitz. For this it is 
enough that Dh be locally Lipschitz. We also need u to be Lipschitz. For this the proofs 
([Ml], [M2], [H]) require only / is a locally Lipschitz twist map whose twist is bounded 
away from zero in a neighborhood of the circle. Again this requires only that Dh be locally 
Lipschitz. 

However, to get a uniform bound in neighborhoods of the lower boundary for the Lip- 
schitz constants of the functions u : S 1 — > (0, 00) that define invariant circles, we need a 
uniform bound on the twist (as defined in Remark 1.1.2). This requires that h\\ and /122 be 
bounded and hence the requirement in (Gl) that Dh is Lipschitz in T> e . 

One can geometrically see the dependence of the Lipschitz constants of the functions u 
by thinking of the twist as defining a cone about the vertical in the tangent bundle at each 
point x. By definition, the image under Df of a vertical vector based at f~ 1 (x) misses this 
cone as does the image of the vertical vector based at f(x) under DQ^ 1 ). This implies that 
vectors in this cone are rotated beyond the vertical in either forward or backward iteration. 
Thus applying the considerations of the previous remark, no tangent to an invariant circle 
can lie in this cone. 

§1.3 Monotone recursion relations. The point of view of Angenent will be used to 
construct examples in §5.1. The basic idea in [A] is to use (1.3) directly to define a recursion 
relation that must be satisfied by sequences of x coordinates of orbits. 

We will state these ideas just in the limited context needed here. Assume that h and hi 
satisfy the properties (Gl) - (G4) above and let 

A(x-i,x ,xi) = -(hi(xo,xi) + fi2(x-i,xo)). (1.5) 

The properties of A that cause it to behave like the monotone recursion relations in [A] 
are: (1) (monotonicity) A is a nondecreasing function of x-\ and x\, (2) (periodicity) 
A(x-\, xq, x\) = A(x_i + 27r,xo + 27r, x\ + 2n), and (3) (coerciveness) which is given by 
the condition in the first part of (G4). The coerciveness condition differs somewhat from 
that adopted in [A], but the proof of Theorem 1.2 below can be constructed with little 
difficulty 

An allowable configuration x is called a solution to the monotone recursion relation given 
by A if A(xi-i,Xi,Xi+i) = for all i. Note that a solution always gives the sequence of x 
coordinates of an orbit of the lift / of the half-cylinder twist map defined by the generating 
function h. The configuration x is a subsolution if A(xj_i, xi, Xi+\) > for all % and a 
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supersolution if A(xj_i, Xi, Xi+i) < for all i. The standard partial order on configurations 
is given by x < y if and only if x% < y% for all i. The main result we need is: 

Theorem 1.2: (Angenent) Let A be defined by (1.5) for a generating function h as- 
sociated with a half-cylinder twist map. If x and y are a subsolution and a supersolution, 
respectively, of the monotone recursion relation generated by A, and further, x < y, then 
there exists a solution z with x < z < y. 

§1.4 The area function of an invariant circle. The next result concerns the con- 
stancy of a certain area that can be associated with an invariant circle. This quantity 
sometimes called "Mather's big H", has a simple geometric interpretation in billiards and 
dual billiards (see the appendix and Theorem 4.2). 

Let f2 be a (not necessarily invariant) circle in C that is the graph of u : S 1 — > (0, oo). 
We will also use u for the function R — > (0, oo) that has Cl C C as its graph. Given the lift 
of a half-cylinder twist map /, we can define two functions associated with the displacement 
of Cl under /. Let a : R — > R be defined as a(x) = 7Ti(f(x,u(x))), or equivalently, a(x) 
is the unique solution to —hi(x,a(x)) = u(x). And let (3 : R — > R be defined by f3(x) = 
7Ti(/(x, <r(x))) where a(x) = sup{7 : (x, 7) G / _1 (fi)}. Equivalently (3{x) = supja;' : 
h,2(x,x') = u(x')} (see Figure 1.3). 




P(x) a(x) 



Figure 1.3: The definitions of a, (3, and the area function. 

From the definitions, a is continuous and (3 < a. Since h\2 < 0, [3 is strictly increasing. 
Note that if 77 > f3(x), then since \im x f^ boo ^ = 00, we have h2(x,r)) > u(rj). Also note that 

fl is /-invariant if and only if a — (3. 

If the generating function associated with / is h, define the area function of f2 as 

f (3(x) 

A(x) := / u(rj) drj - h(x,(3{x)). 

Thus A(x) measures the signed area between the Q and the image of the vertical arc above 
x (see Figure 1.3 again). Mather proves in [M3] that the area function of an invariant circle 
is always constant. 

Proposition 1.3: (Mather) Given a half-cylinder twist map f , the graph Q of a con- 
tinuous function u : S 1 — > (0, 00) is an f -invariant circle if and only if its area function is 
constant. 
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Proof: First assume that Q is /-invariant. By Theorem 1.1(b), u is a Lipschitz function 
and since / is locally Lipschitz, a = (3 is Lipschitz and so its derivative exists almost 
everywhere. Differentiating A(x) using its definition and then applying (1.2) one gets that 
dA/dx = almost everywhere, and since A is continuous, it is constant. 

Now conversely, assume that A is constant and proceed by contradiction. If a ^ (3, then 
since a is continuous and (3 < a, there exists x and x with x < x and at{rj) > (3(x) for all 
7] G (x,x). Thus for these rj, hi(rj, (3(xj) > hi(r],a(r])) = —u(rj). Note also that above we 
showed that h2(x,r)) > u{rj) for all i] > j3{x). Now by assumption, A(x) — A{x) = 0, and so 



rx _ _ rpyx) 

— j u(rj) dr] = —h(x,P(x)) + h{x ) f3{x)) — j u{rj) dr] 
= ~h(x,p(x)) + h{x,P(x)) - u(r]) dr] 

Jf3(x) 



- h(x,P(x)) + h(x,P(x)) 

P{x) rX 

(h 2 (x, rj) - u{rj)) dr] + / hi(x, f3{x)) dr] 

(3{x) Jx 



> 

I X 



f x 

/ u{rj) dr] 

Jx 



where in the last inequality we have used the fact that (3 is strictly increasing. □ 

§2 Convex curves and envelope coordinates. 

In this section we review some standard material from the theory of convex curves in 
the plane. For more details the reader is referred to [E], [Gl], and [G2]. 

§2.1 The height function and the radius of curvature. If V is a convex closed curve 
in the plane, let U(T) denote the unbounded component of M 2 — V. It will be convenient to 
assume that the origin is not contained in U (T) and that Y is oriented in a counterclockwise 
direction. Let ug = (cos(0), sin(0)) and n# = (- sin(0), cos(0)). For each 9 e S 1 = R/(2ttZ), 
£.0 is the supporting line to T that is parallel to n# and near its point of intersection T is 
oriented in the same direction as n#. It will also be convenient to assume that Co intersects 
T in an extremal point. The height function of T, p(0), is the distance from the origin 
to Cq (see Figure 2.1(a)). Equivalently, p : S 1 — > K. is given by p(9) = G(ug) where 
G(z) = maxj^ • w : w G T}. 

If Cq intersects T in an extremal point, let a(8) be this intersection. If Cg intersects 
T in a non-trivial segment, let a(8) be the extreme endpoint of the intersection in the 
counterclockwise direction. The arc length along T from a(0) to a{9) is denoted by s(0). 
Note that s is nondecreasing, continuous from the right and 

,0 

a{9) — a(0) = / n ri ds(r)) 

J 

as a Riemann-Stieltjes integral. The points of discontinuity of s are exactly the countable 
set of for which Cg intersects T in a non-trivial segment. 
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Figure 2.1: (a) The height function; (b) The geometry 
of equation (2.1). 



A simple geometric argument (see Figure 2.1(b)) shows that 



p(9) = -n e ■ (a(0) - a(0)) + q(0) sin(0) +p(0) cos(#) 




(2.1) 



In Figure 2.1(b) the three terms in this expression are labeled A, B and C, respectively. The 
quantity q(0) is as pictured, explicitly, q(0) is the signed distance from a(0) to p(0) uo (the 
q(0) shown in Figure 2.1(b) is positive). Note that the last two terms in (2.1) just reflect 
the position of the origin. 

It follows from (2.1) that p is differentiable at all but a countable number of points where 
p' has a simple jump discontinuity and 



At points where p' has a jump discontinuity (those 9 for which Cq intersects T in a segment) 
this formula may be interpreted in terms of right and left hand limits. In particular, g(0) = 
p'(0) and p' has a geometric interpretation as the signed distance from a{9) to p{9) ug (see 
Figure 2.1(a)). We shall be primarily interested in the case when V is strictly convex which 
corresponds to the condition that s(9) is continuous and p is C 1 . 

If the arclength function s(9) is differentiable, p{9) := ds/d9 is the radius of curvature 
of T, i.e. the radius of the osculating circle to V at the point a(9). Using (2.1) and (2.2), 



If the s{9) is not differentiable, this equation may be interpreted in terms of distributions. For 
example, if Y is a polygon, p is the sum of weighted Dirac delta functions. In any case, since 
s is nondecreasing, p exists almost everywhere. One interpretation of (2.3) is that the height 
function is given by a forced harmonic oscillator with forcing function given by the radius 




(2.2) 



P + P = P > 0. 



(2.3) 
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of curvature (in which case the last two terms in (2.1) reflect the initial conditions). Very 
informally, a closed curve is a point that has been forced outward by the radius of curvature. 
This observation will be crucial in connecting dual billiards with impact oscillators in §7. 

Remark 2.1.1: For later reference we give a characterization due to J. Green of the 
type of periodic functions that can arise as height functions for a convex closed curve in the 
plane ([Gl], [G2]). The characterization involves a generalized convexity. Let S(x; A, B) = 
A cos x + Bsinx. A function g : M. — > 1R is said to be sub-sine if for each x\ < X2 < x\ + 7r, 
S(x; Aq, Bq) > g(x) for all x G [xi,X2], where Aq and Bq are the unique values for which 
S(x\; Aq, Bq) = g(xi) and S(x2] Aq, Bo) = g{x2}- Green shows that a 27r-periodic function 
p is the height function of some convex closed curve if and only if p is sub-sine. This 
happens if and only if p can be written as in (2.1) for some nondecreasing function s with 
/ exp(i#) ds(8) = 0. Further, a twice differentiable p is sub-sine if and only if p" + p > 0. 

In a similar vein note that if s : [0, 2tt) — > M. is any bounded, nondecreasing function 
with J exp(i9) ds(9) = we can construct a T using (2.1) and (2.2). The resulting curve 
will have parameterization a{9) = Rg(p(9),p'(9—)) where Rg is rotation in the plane by the 
angle 9 and p'(9—) is the derivative at 9 from the left. The curve T will be the convex hull 
of a(S 1 ). One can also use (2.1) to identify convex curves whose perimeter has length one 
with regular Borel probability measures ds on S* 1 whose first Fourier coefficient vanishes. 
Atoms in the measures give flat spots (intervals in T); strictly convex curves correspond to 
measures without atoms. 

§2.2 Envelope coordinates. The next step is to use T to give coordinates to the open 
(topological) annulus U(T). For a point given in Euclidian coordinates z G U(T), let 9(z) be 
the unique angle 9 with z G Cg and a(9) — z is parallel to n#. A second coordinate for z is 
given by the distance from z to a(9) along Lqi z \ and is denoted £(z) (see Figure 2.2). Note 
that the area elements are related by dz = £ d£ d9 = d'j d9 where 7 = £ 2 /2. Thus the map 
z I— > (9, 7) is area preserving. The coordinates (9, 7) will be called envelope coordinates. If 
T has corners, i.e. its radius of curvature function p{9) vanishes on a nontrivial interval, the 
map z I— > (0,7) cannot be extended to a single valued function on Y. 



§3 Dual billiards yields a half-cylinder twist map 

A geometric description of the dual billiard map arising from a convex curve Y was 
given in the introduction. The purpose of this section is to show that under appropriate 




z = (6,Y) 



Figure 2.2: Envelope coordinates. 



13 



assumptions on T the dual billiards map in envelope coordinates is a half-cylinder twist 
map. It is geometrically clear that a dual billiards map preserves Euclidian area. The 
corresponding map in envelope coordinates is area-preserving because the change to envelope 
coordinates preserves area. The other main property that characterizes twist maps is that 
the image of a vertical ray is a graph over an arc in the circle. A vertical ray above the point 
(6q, 0) in envelope coordinates corresponds to the negative ray in Cg beginning at the point 
of tangency with T. The image of this ray under the dual billiards map is the positive ray 
in Cg Q that begins at the same point. Thus the condition that makes / a twist map is that 
this ray should hit each other Cg (9o < 9 < 9q + 7r) in exactly one point. This follows from 
the convexity of T. 

The dual billiards map on U(T) is denoted §y> or just an d the corresponding map in 
envelope coordinates is denoted /r or /, and its lift to C as /. 

§3.1 The recursion relation from geometry. Let T be a convex curve in the plane. It 
has associated functions a, p, etc. as defined in §2.1. Given two angles 9 < 9' , let 5 = 9' — 9 
and = a{9') — a{9). If R = R(9, 9'), L = L{9, 9'), u\, and u>2 are as pictured in Figure 3.1, 
then the law of sines yields 

sin(7r - 5) _ sin(c^i) _ sin (u 2 ) ( „ , , 

|0l ~~ I—'— R~- [6A) 




Figure 3.1: Geometric basis of the recursion relation. 



The wedge (or cross product) between two vectors in the plane is w A v = W1V2 — v 2 wi = 
w\\v\ sin(c<j) where uo is the angle from w to v. Using this notation and (3.1) we have 



r(f) 9') l^l sin (^l) ng A f e n v ds{r]) _ J e sm(r] - 9) ds(r]) 

{ ' ' sin(7r - (9' - 9)) sin(0' - 9) sin(0' -9) ' 1 ' ' 



and a similar calculation yields 



me)- — 8in(6)/ _ e) — • (3-3) 
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Thus using the geometric definition of the dual billiard map a triple of points (9, 9', 9") are 
the angles of an orbit of dual billiards if and only if < 9' — 9 < n, < 9" — 9' < ir, and 

L(9,9') = R(9',9")>0. (3.4) 



§3.2 The generating function and the boundary maps. Next we define the functions 
h and bi that are the generating function and boundary maps for /p. If we let 



h(x, x) = - (L(x, y)Y dy = - (R(y, x')) z dy (3.5) 

Z Jx Z Jx 

then hi = —R 2 /2, h<i = L 2 /2, and the condition that h\{x' ,x") + h,2(x,x') = is equivalent 
to (3.4). 

To define the boundary maps, let bo(x) = inf{z : z > x and p(z) > 0}. Perhaps put 
more simply, bo(x) = x if x is not in an interval in which p identically vanishes. If x is 
in such an interval, bo(x) is the right endpoint of that interval. The function bo is clearly 
continuous from the right, nondecreasing, satisfies bo(x + 2n) = bo(x) + 2tt and has rotation 
number equal to zero. Let boo(x) = x + n. 

Let us assume now that T is a convex curve whose arclength function s(0) is Lipschitz. In 
this case the radius of curvature p exists almost everywhere, is bounded and satisfies p > 0. 
The derivatives at points where p exists are: 

Ri(x, x) = —p(x) + cot(x' — x)R(x, x') 
, L(x,x') 

R2{X,X) = —— 

sm{x' - x) 

T ( , RM (3 - 6) 

Li(x,X) = — r 

sm(r — x) 
L>2(x, x) = p(x') — cot(x' — x)L(x, x'). 



It is now a simple matter to check that h and the bi satisfy the necessary properties to define 
a half-cylinder twist map. 

Theorem 3.1: If T is a convex curve whose arclength function s(6) is Lipschitz, then 
its dual billiard map /r in envelope coordinates is a half-cylinder twist map with rotation set 
equal to [0, 1/2]. 

Remark 3.2.1: The function h(9,9') has the geometric interpretation as the area 
bounded by T and the supporting lines Co and £g> ([Ms4]) (see Figure 3.2). This area 
is the same as the geometric interpretation of the generating function given in §1.1 because 
the change from Euclidian to envelope coordinates preserves area, The variational formula- 
tion also has a nice geometric interpretation in dual billiards. For example, to find a period 
three orbit of the dual billiards map, one finds the circumscribed triangle of smallest area 
(see Figure 3.3, left). Theorem 1.1(a) says that we can find a monotone periodic orbit of 
all rotation numbers. This corresponds to finding circumscribed polygons of each "rotation 
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Figure 3.2: The geometric interpretation of the action. 




Figure 3.3: Using variational methods to find mono- 
tone periodic orbits with rotation number 1/3 (left) and 
2/5 (right). 



type" (see Figure 3.3, right). Further, there are the irrational analogs of these sets for each 
irrational in the rotation set. 

Remark 3.2.2: If the curve T is strictly convex, then its arclength function s(9) is 
continuous, and the half-cylinder twist map /$ on C is also. By inspecting (3.6) we see that 
this is sufficient to make h\2 continuous, even though tin, etc. may not exist. If s is C r , 
then so is /p. 

If s is Lipschitz, Dh will be also, and so is locally Lipschitz. However, as \x — x'\ — > 0, 
hi2 — > 0, and so (1.2) yields that / is never Lipschitz, no matter how smooth s is. On 
the other hand, / has bounded twist near the bottom boundary because h\\ and hi2 are 
bounded there. In fact, as \x — x'\ — * 0, h\\ — > 0, so the slope of the image of a vertical 
tangent goes to zero near the bottom edge. 

Note that $r always extends continuously to T by making it the identity there, but /r 
will not have a continuous extension to S 1 x {0} if T has corners (nontrivial intervals where 
p vanishes). This is why discontinuous boundary maps &o had to be allowed in the definition 
of half-cylinder twist map. Note also that as x' — x — > n (i.e. as 7 — > 00), /122 00 so the 
image of vertical tangent is near vertical and the twist of /p goes to zero. 

§4 Invariant circles for dual billiards. 

As noted in the introduction, the question of the existence of invariant circles is central 
in twist map theory. An invariant circle for a dual billiards map is geometrically connected 
to the given convex curve via an equal area construction. Recall that in this paper the phrase 
"invariant circle" always means a homotopically nontrivial circle. 
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§4.1 Existence of invariant circles. The first theorem gives sufficient conditions on the 
curve T so that its dual billiards map has invariant circles near to T and near infinity. The 
Hausdorff distance between two compact sets, X and Y, is denoted d(X, Y). 

Theorem 4.1: (R. Douady) Let Y be a convex curve in the plane with a radius of 
curvature function p{6) that is C r . 

(a) If r > A, then for each N > the set of <3>r '-invariant circles Q with d(Q, T) > N 
has positive Lebesgue measure. 

(b) If r > 5 and p > 0, then for each e > the set of $p-mwarmni circles Q with 
d(Q, r) < e has positive Lebesgue measure. 

The basic idea of the proof in [D] (see also [Ms2]) is to extend the half-cylinder twist 
map / associated with <3>r to the boundaries of C* and then use a normal form result based 
on Herman's version of the curve translation theorem ([H]). In the compactification of the 
half-cylinder, / does not preserve a finite measure near infinity. However, it does satisfy the 
circle intersection property, so the curve translation theorem can be used to obtain invariant 
circles. The condition p > in (b) is needed to extend / smoothly tolx {0} (c/. Theorem 
4.3(a) below). 

§4.2 Invariant circles and area envelopes. The next result concerns the inverse problem 
for invariant circles, i.e. given a closed curve Ti can you find a convex curve To so that Ti 
is an invariant circle for $r - The solution involves a geometric construction called the area 
envelope. 

Given a simple closed curve T\ in the plane and a number a less than half the area 
enclosed by T\, for each angle we can find a unique oriented line Cg in the direction ug 
so that the area to the right of the line and inside Ti is equal to a. Denote the envelope 
of these lines as AE(Ti,a) (see Figure 4.1). If AE(Ti,a) is convex, it is an easy geometric 
exercise to show that the locus of the midpoints of the chords of the lines Cg in Ti is the 
area envelope." The next result is a consequence of this geometric fact. It appeared first in 



[D]. 




Figure 4.1: An area envelope. 



" The author learned this from R. Ticciati who "discovered" it during a High School project 



in 1967 
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Theorem 4.2: (R. Douady) If Y is a convex curve with a Lipschitz arclength function, 
then a closed curve Y\ C U (Y) is invariant under the dual billiards map $r if an d only if 
Y = AE(Ti, a) for some a. 

Recall that the change to envelope coordinates is area-preserving. Thus if one considers 
the half-cylinder twist map corresponding to $r> the quantity A(x) (defined in §1.4) asso- 
ciated to the invariant circle is exactly the appropriate area a for which T = AE(Y\,a). 
Thus Theorem 4.2 is a consequence of the more general fact given in Proposition 1.3. The 
appendix contains a remark on an analogous result for inner billiards. 

Even if Ti is convex, its area envelope need not be (eg. an area envelope inside an 
equilateral triangle). With the appropriate interpretation of chord, the area envelope will 
still be the locus of the midpoints. Area envelopes are discussed in more detail in [FT] and 
[GK]. 

§4.3 The non-existence of invariant circles. The next result gives sufficient conditions 
on T that imply the non-existence of invariant circles near T. These results are the analogs 
for dual billiards of results for billiards due to Mather and Hubacher. Conditions on the 
curvature in billiards play the role of the conditions on the radius of curvature in dual 
billiards. 

The proof of part (a) goes exactly like that of Mather in [Ml], and indeed, Theorem 1.1(c) 
is taken from that paper. The proof of (b) is unlike that of [Hr] in that we use estimates 
involving the recursion relation of the twist map, rather than geometric arguments. Part (a) 
is also contained in [GK]. That paper also contains estimates based on the geometry of T for 
the sizes of regions in U (Y) that do not contain invariant circles. 

Theorem 4.3: Let Y be a convex curve in the plane with Lipschitz arclength function 
s{9), and J <E S 1 is an open interval. If 

(a) the arclength function s(9) is C 1 on J, p> at the endpoints of J, and for some 
6eJ, p(9) = 0, 

or 

(b) there is a 9 G J so that s is C 1 on J— 9, p > on J , and p has a jump discontinuity 
at 9, 

then there exists a neighborhood ofY that contains no homotopically nontrivial $r-wyaha??i 
circles. 

Proof: The proof will use the half-cylinder twist map / associated with Y that has lift 
/. Let h, bo and V be the generating function, lower boundary map and the domain of h, 
respectively, and x is a lift of 9. 

If Q C C is the lift of a homotopically nontrivial /-invariant circle, then by Theorem 
1.1(b), Q is the graph of a Lipschitz function u : R — > (0, oo). Using the derivative formulas 
(3.6) and Theorem 1.1(d), for almost every xo, 

2p(xo) > cot(xi — xo)R(xo, x\) + cot(xo — x-i)L(x-i,xo) (4.1) 

where (a;;, 7;) = f(x , u(x )). 
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Assume first that the situation in (a) happens. If there were a ^-invariant circles arbi- 
trarily close to to T, then by Theorem 1.1(b), we could find the lift of an invariant circle Q ar- 
bitrarily close to Ix {0} in the Hausdorff topology, in which case max{|a;i — xq\, \xq — x-i\} < 
n/2 and {x-i, xq, x{\ C J for all xq sufficiently close to x. In this case the right hand side of 
(4.1) is positive, continuous and bounded away from zero while the left hand side is a con- 
tinuous function that vanishes at xq = x. Thus the equation cannot hold almost everywhere 
for xq in a neighborhood of x, a contradiction. 

Now assume the situation given in (b). Let lim^^i p(x) := p ± and assume without loss 
of generality that p + > p~ and \J\ < n/2. The principle estimate needed is the following: 
Ve > 0, 35\ > 0,ViV > , 36 2 = 62(61, N) so that \x-i - x \ < 6\, 6\/N < \xq-x\\ < 5\ and 
< x — xq < 62 implies 

max{|Z(a;_i,a;o) ~ P~\, \R(xq,x{) - p + \, \p(x ) - p~\} < e 

where L(x-\,xq) = cot(xo — x-i)L(x-i, xq) and R(xq,x\) = cot(xi — xq)R(xq,x\). 

To prove this first note that if (x-i, xq, x\) is a segment of an allowable configuration 
for h, {x-i, xq, xi} G J, and x > xq, then 

. - , N 1, , x f Xl , sin(xi — 77) tan(xi — xq) , . 

\R(x-i, xq) - p^ = cot(xi - xq) (— ; -p(r]) s-P ) dr l\ 

Jx sm(x\ - xq) (xi - xq) 

1 r x i 
< 7 r / \p(v) ~ P I di] 



and 



< 7 r((x - x )sup{\ p(rj) -p + \:r]e [x ,x)} 

(xi - xq) 

+ (xi - x)sup{\p(r]) - p + \:i] E (x,xi]}) 

1 f x o _ 

\L(x- 1 ,xq)-p \ <- -/ \(p(rj)-p )\dr] 

(Xq — X--I 



1 r x o 

Xq — X-l) Jx-i 

< sup{\p(r]) - p~\ : rj e [x_i,x ]}. 



Now given e > 0, since p is continuous and bounded on J — x there is an 6\ > so that 
< x\ — x < 6\ implies 

sup{|p(r?) - p + \ : 77 e (x,X!}} < e/2 
and in addition, < x — X-\ < 26\ implies 

sup{|p(r/) - p~\ : 77 G (x-i,Xq)} < e. 
Now given N, pick 62 < 6\ so that 

sup{|p(?7) — p \ '■ f] E [xq,x)} < e/2. 



61 

With these choices, note that x — x-i = xq — x-\ + x — xq < 61 + 62 < 26\, and thus principle 
estimate follows. 
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Continuing the proof of the lemma, pick e > so that p + — p > 4e and let Si be as in 
the principle estimate. If there are /-invariant circles arbitrarily near S 1 x {0}, then there is 
a lift of an invariant circle for which < x% — x%-i < Si for all X{ G J where X{ = iri(f n (z)) 
with z an element of the lift of the invariant circle. Fix one such circle Q. Since it contains 
no fixed points by Theorem 1.1(b) there is an N with Si/N < Xi — Xi-i < Si for all pairs 
on Q. Using this N, find S2 as in the principle estimate. Then the principle estimate along 
with (4.1) implies that for almost all xq, in particular for some xq with < x — xq < 62, 
> L(x-i,xq) + R(xq, xi) + 2p(xo) > p + — p~ — 4e > 0, a contradiction. □ 

Remarks: 

4.3.1: Neither the condition of part (a) nor (b) excludes the existence of invariant circles 
outside neighborhoods of T. Indeed, the area envelope construction of §4.2 allows one to 
construct an invariant circle for dual billiards on a convex curve whose radius of curvature 
function has any prescribed local behavior. In the case of billiards, the analog of part (a) (a 
point where the curvature vanishes) excludes all invariant circles, while the analog of part 
(b) does not. 

4.3.2: The Birkhoff-Mather stability theorem relates the non-existence of invariant cir- 
cles with the existence of orbits with certain limit behavior. The theorem states that if Hi 
and Q2 are invariant circles for a twist map, then there are no other invariant circles in the 
annulus bounded by Qi and Q2 if and only if there is a point z whose a- and u;-limit sets are 
contained in Qi and Q2, respectively. Mather's result is proven in [M4] for a specific class of 
twist maps, but the arguments are quite robust and almost certainly apply to half-cylinder 
twist maps. 

Thus, for example, when there are no $p- invariant circles near T there is an orbit that 
converges to T in forward time. This orbit is distinguished from the crash orbits of the next 
section by the fact that it converges to all of T, not a single point on T. There will also be 
an orbit that converges to T under backward iteration, so it may be viewed as "emerging" 
from T. 

The question of stability depends on the behavior at infinity. In the Birkhoff-Mather 
theorem, one of the invariant circles can be located at infinity. Thus if there are no invariant 
circles near infinity, there is an orbit that escapes, i. e. goes to infinity under forward iteration. 
As noted in the introduction one of the most important outstanding problems in the theory 
of dual billiards is the existence of a convex curve whose dual billiards map has such an 
orbit. 

4.3.3: Remark 1.2.1 contains the main geometric idea underlying the proof of Theorem 
4.3. The conditions given in (a) or (b) ensure that a vertical tangent tumbles past the 
vertical in two iterates. To obtain a stronger result on the non-existence of invariant circles, 
one would need to use more iterates. However, for the half-cylinder twist maps coming from 
dual billiards, the twist goes to zero at infinity. It thus requires more and more iterates 
to tumble a tangent as one approaches infinity. This partially explains the difficulty in 
obtaining a T that yields no invariant circles near infinity. As noted in Remark 4.3.1, no 
local condition on T can eliminate these invariant circles. 
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§5 Crash orbits 

A crash orbit for dual billiards is an orbit that converges under forward iteration to a 
point on the convex curve T. The first theorem states that there exist convex curves with 
differentiable radius of curvature functions p for which there are crash orbits. However, the 
second theorem says that if p is positive and has bounded derivative, then there do not exist 
crash orbits. These are analogs of billiards results due to Halpern ([Ha]). The method of 
proof here is different, although the analog of Hapern's construction would also suffice. 

We give a fairly explicit p and then use the criterion of Angenent from Theorem 1.2 to 
get the existence of the crash solutions. The simplest version of the example has a radius 
of curvature function near zero that looks roughly like 1 + (l/2)x 2 sin(l/x 2 ). Note that this 
is a standard example of a function whose derivative exist everywhere, but the derivative is 
not bounded. 

The existence of the crash orbits implies that the dual billiards map has no invariant 
circles near T. However, as in Remark 4.3.1, one can construct examples with crash orbits 
and invariant circles far away using the area envelope construction. 

§5.1 Examples with crash orbits. 

Theorem 5.1: There exists a convex curve T with a differentiable radius curvature 
function such that the corresponding dual billiards map has crash orbits. 

Proof: We use the Angenent criterion of Theorem 1.2 with a modified form of the 
recursion relation. Define A(x_i, xo, x\) = R(xq,xi) — L(x-i,xo) where R and L are as in 
(3.2) and (3.3). It is clear that solutions, subsolutions, etc. of this recursion relation will be 
the same as those given by the A defined in (1.5) with h given by (3.5). 

We will construct a subsolution c with c„/0, and a supersolution w with w n / w < oo, 
so that c < w. By Theorem 1.2, there exists a solution x with c < x < w, and since any 
solution is of necessity monotone increasing, x n converges as n — > oo. Now by (1.1), there 
exists some (#0,70) w hh %i = ^i(P(xo, 7o))- Then using (3.2) or (3.3), 7 n — > 0, and so 
/*(#o,7o) converges to a point on R x {0}, thus the corresponding orbit in U(T) converges 
to a point of T under the dual billiards map. 

Given a sequence a n (to be specified later) with n G N such that a n /* and a n G [— .1, 0), 
let S n = (a n +i — a n )/2 and b n = a n + 5 n . Pick a < c < .001, k > 1, and for n — 1, 2, . . . , 
define p : R -> R by 



p{x) = < 



1 — ca^ for x G [b n -i, a n ) 
1 + ca^ for x G [a n , b n ) 
1 for x G [0, 2tt - b ) 



and extend p so that it is 27r-periodic. 
Now let 

r(x) = 

and 

s(x) 



cos(x) — cos(2a;) 
sin(2a;) 

1 — cos (a;) 
sin(2x) 
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Note that at points where the denominator is zero, r and s can be continued to real analytic 
functions with s(x) +r(x) = tan(x). 

If r n = r(S n ) and s n = s(S n ), a simple calculation yields that 

R n := R(a n , a n+1 ) = (1 + ca k n )r n + (1 - ca n+1 )s n 
L n -i := L(a n -i,a n ) = (1 + ca^-^Sn-i + (1 - ca*)r„_i. 

Our goal is to show that R n — L n -i > for sufficiently large n. Assuming 5 n -i — 5 n >0, 
this is equivalent to 

c ^ a n r n Q"n-\-iSn a n —\S n —i + a n T n —\^ 

_ s n—l ~\~ r n—l s n r n tr t\ 

> x x ( 5-1 ) 

_ tan(^ n _i) - tan(£ n ) 

0~n-l - 0~n 

which as n — ■» oo goes to tan'(O) = 1. 

Now specify a n = —n b for some b < 0. In this case, 

-bn b - 1 >5 n > -bin + l) 6 " 1 
b{b - l)(n + l) b - 2 < <5 n _i - S n < 6(6 - l)(ra - 1) 6 ~ 2 . 

Using Taylor's theorem, for sufficiently small positive 5, r(<5) > 35/4 and —s(S) > —5/2. 
Thus the left hand side of (5.1) is larger than 

\n kh {{n + l) b - 1 + n h - 1 ) -\{n+ l)^" 1 - j(w - l)^ 1 ) 6 " 1 
C( (l-6)(n-l) 6 - 2 } ' 

which goes to infinity if kb + 1 > 0. Thus for these choices of fc and 6, i? n > L n -i for n > N 
for some sufficiently large N. 

Now pick p that is C°° except at and is close enough to p to ensure that a is still a 
subsolution (for n > N) for the recursion relation defined using p. 

To construct the subsolution c, let 7jv = ^2( a iV-i) a iv) and for i < 0, c^r+i = P(o>n,1n) 
and for i > 0, c^+i = Gtjv+i- By construction, A(cj_i, Cj, Q+i) = 0ifi<A^ — 1 and 
A(q_i, q, Cj+i) > if i > TV, and thus c is a subsolution. 

To construct a supersolution w with w > c begin by letting V{ = 7ri(/*(a/v, 7iv)) where 
(ajv, 7jv) is defined in the previous paragraph. Now if the sequence t> j is bounded above it is 
a crash orbit and we are done, otherwise let iq = inf {i : V{ > 0} and w^^j = Vi ^j for j > 0. 
Define w^+j for j > inductively as follows. Assuming w n -\ and w n have been obtained, 
pick w n < w n+ i < wn + T,i=o~ N 2»T3 so ^at < R(w n , w n +\) < L(w n -i, w n ), where L and 
R are defined by (3.2) and (3.3) using p for the radius of curvature. This is possible because 
R(x, x') — > as \x' — x\ — > 0. 

Now note that A(tUfc_i, u>fc, ttffc+i) is less than zero for A; > iV and equal to zero for 
k < N, and thus w is a subsolution. By construction, c < w and limw n < wn + 1/4. □ 
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§5.2 Nonexistence of crash solutions. 

Theorem 5.2: If T is a convex curve in the plane with a radius of curvature function 
p > that is differentiate and the derivative is bounded, then the dual billiards map has no 
crash orbits. 

Proof: Assume to the contrary that the dual billiards map has a crash orbit. Lift 
the corresponding half-cylinder twist map to the universal cover and say the crash orbit is 
(x n ,ln), where we assume without loss of generality that x n S 0. 

Using the recursion relation given by (3.4) and the mean value theorem for integrals we 
can find i] n ^i G (x n -i,x n ) and r n G (x n , x n+ \) with 

tan(5 n _i/2) = p(r n )tan(5 n /2) 

where 5 n = x n +i — x n , and thus 

(p(Vn-i) ~ p(i~n)) tan(5 n _i/2) = p(r n )(tan(5 ri /2) - tan(5 n _i/2). 

Now by construction, r n — i] n -i < 8 n + 5 n -i and if 5 > is small, 5 < tan(5) < 25. 
Using the mean value theorem for derivatives, there is a a n G [%-i, r n ] with 



\p'M 



> 



| PK) -p(ri n -i) 
T n ~ Vn-1 

p(r n )(tan(5 n /2) - tan(^ w _i/2)) 

(r n -%-i)tan(5 n _i/2) 
| p(r n )(^ t /2-^_ 1 ) | 
(o~n + 0~n-l) 0~n-l 



2 5 n _i + 5 n 

Now since p is continuous, as n — > oo, p(t„) — > p(0) > 0; since X)^n converges, <5 n — > 
and limsup |5 n /5 n _i| < 1. Thus the last expression in (5.2) diverges, contradicting the 
boundedness of p' . □ 

§6 The normalized action of the derivative 

In this section we make some geometric observations about the action of the derivative 
of the dual billiards map on tangent directions. Since our aims are primarily descriptive, we 
will adopt a more informal tone than in previous sections. 

§6.1 Computation of the action of the derivative. Let (3(t) = (/9i(t), foit)) be a 
parameterization of Y by arc length, and let us assume initially that (5 is twice differentiable. 

Fix a point z G u(Y) and say its envelope coordinates are (9, 7). Let t and / be such that 
f3{t) = a{9) and t = \^2~y. Thus j3(t) is the point on T used for the dual billiards shot, and 
£ the distance from z to j3{t). From inspection of Figure 6.1 we have 

2(3{t) = z + $(z) (6.1) 
z = (3(t)-£(3'(t). (6.2) 
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*(z) 




Figure 6.1: Some geometry of dual billiards. 

Adopting the classical (and often confusing) notation of using the same symbol for a 
coordinate and a coordinate change, z = (x,y) = (x(t,£),y(t,£)), and subscripts to denote 
differentiation, (6.2) yields 



M :-. 



x t xi 

yt yi 



P[-m Pi 
$2-02 & 



Note that det(M) = £((3' A (3") = tn(t), where n(t) is the curvature at the point (3{t) (the 
curvature is the reciprocal of the radius of curvature p). Thus inverting M we get 



(t x 


ty) 


nt ( 


\£ x 







P'2 



-P'l 



(6.4) 



Differentiating (6.1) yields 

D$ + Id = 2 



f3[ 



k£\{^ 2 



and thus acting on a vector u, 



t X ty 

£ X ty 

-{P'lf 



(D$ + Id)[u] = -(u A (3') (3'. 



(6.5) 



We now need some elementary geometry of the wedge product. Let w be a unit vector, 
u and v be vectors, and a, b, c be positive scalars so that bw = av + cu. If A represents the 
area of the triangle formed by the vectors then, 

2A = cu A bw = bw A av = cu A av 



and so 



2A 



b 2 (w A v)(w A u) 
v A u 



(6.6) 
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Now wedging (6.5) by £>$[it] on the left yields 



W[u]Au = -(uA P')(D$[u] A f3') 



(6.7) 



Thus if B is the area of the shaded region in in Figure 6.2, i. e. B is the area of the triangle 
bounded by the line connecting z to $>{z), the line in the direction u passing through the 
point z, and the line in the direction of -D$[u], then using (6.6) and (6.7), 



= (2f) 2 (^AJ$[u])(^Au) = ^ 
2 D$[u] A u 

Equation (6.8) contains the same information as Proposition 2.2 in [GK]. 




Figure 6.2: Geometry of the derivative. 



(6.8) 



§6.2 Geometric computation of the derivative. To geometrically compute the action 
of the derivative of the dual billiards map on a tangent direction, first fix a point z and 
its image $(z) as in Figure 6.2. Now draw a line parallel to that from z to $(z) that is a 
distance k£ 2 from that line, where k is the curvature at the point of tangency. Call this new 
line the bounce line. Given a unit tangent vector v based at z, draw a line in the direction 
of v until it hits the bounce line. The line from this point of intersection to the point &(z) 
will give the direction of the vector D<E»[v]. This is because the shaded region in Figure 6.2 
has area equal to k£ 3 as required by (6.8). 

If T has a point with p{6) = 0, by taking limits one sees that the action on tangent 
directions is a bounce off infinity. Thus when the dual billiards shot uses a corner of T, the 
action of _D$ is rigid rotation by it. 

Remark 6.2.1: One motivation for seeking a geometric understanding the action of the 
derivative is to assist in the finding of strictly invariant cone fields and thus (perhaps) obtain 
positive Lyapunov exponents as in [W]. This theory does not seem to be well understood in 
the case of infinite area as is the case in dual billiards. One could also search for these cone 
fields in bounded regions, eg. between a given curve and a second curve inside constructed 
using the area envelope construction. The existence of positive exponents is particularly 
interesting in light of the physical interpretation of dual billiards given in §7. 

§6.3 Comparison with a hyperbola. Having understood geometrically the action of the 
derivative we are now in a position to interpret the tumbling tangent criterion of Remark 
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1.2.1. Again, pick a point z. Figure 6.3(a) shows the action of the second iterate of a 
dual billiards map on a tangent direction based at &~ 1 (z). There are three supporting lines 
involved in the iteration, call them C-i,Cq, and C\. The half-ray in C-\ that begins at Y 
and points downward corresponds to a vertical line in envelope coordinates. Thus the initial 
vector chosen corresponds to the initial vertical tangent in Figure 1.2. Its first iterate under 
yields a parallel vector pointed in the opposite direction. For the next iterate we use 
the procedure of the previous subsection. Draw a bounce line at height k£ 2 and compute 
the next image. The condition that this vector has tumbled beyond the vertical is that the 
image is on the side of L\ nearest V. In the figure the initial vector has tumbled by the third 
iterate. 




Figure 6.3: (a) Tumbling tangents in dual billiards; 
(b) The comparison hyperbola. 



It is geometrically clear now that when the curvature is infinite at the point of tangency 
of Cq (i.e. the radius of curvature vanishes) one always gets tumbling tangents for points z 
near T, and thus no invariant circles near T. It is also clear why the same argument does 
not work for z that are far away from T. 

To get a more refined criterion, construct a hyperbola as the area envelope inside the 
pair of lines £_i and C\ using an area equal to that of the triangle bounded by the C-i, jCq, 
and C\ (see Figure 6.3(b)). This hyperbola will be tangent to Cq at the same point that Y 
is. Now if \1> is the dual billiards map of this hyperbola, the image of tangent vector based 
at z that points in the direction of £_i will be a vector based at $(z) that points in the 
direction of C\. Thus by §6.1, the point £_i n C\ is a distance k p £ 2 from Cq where k p is the 
curvature of the hyperbola at the point of tangency. Recalling the geometric condition for 
a tumbling tangent for $ one sees that the tangent tumbles if and only if the bounce line 
for T is higher than the bounce line for the hyperbola, i.e. if ky > k p . A simple calculation 
shows that this condition is identical to (4.1). 
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§7 Dual billiards and impact oscillators 

In this section we show that the dynamics of a dual billiards map associated with a 
convex curve is identical to the dynamics of a certain impact oscillator. As in the previous 
section, our purpose is primarily descriptive so we maintain a somewhat informal tone. 

§7.1 The impact oscillator. The first ingredient in the impact oscillator is a periodically 
moving wall. Its position at time t is given by pit), and we require that p(t + 2n) = pit), 
p is continuous, and p is sub-sine as defined in Remark 2.1.1. As noted there, if p is twice 
differentiable this happens if and only if 

p + p:= p > 0. (7.1) 

Alternatively, we can start with a periodic p with / pit) exp(ii) dt = 0, pick initial conditions 
for p{0) and p{0), and then solve (7.1) for p. 

The second ingredient in the impact oscillator is a particle whose position is given by 
x(t). Except at collision this particle is a simple harmonic oscillator with period 2n, i.e. it 
satisfies 

x + x = 0. (7.2) 

The collisions of the particle with the wall are assumed to be perfectly elastic and the 
wall has infinite mass; so if t c is a time when x(t c ) = p(t c ), then 

-(xbeforeiQ ~ p(*c)) = iafteAQ ~ p(t c ). (7.3) 

A solution of the impact oscillator is a function x(t) > p(t) that satisfies (7.2) except 
when x{t) = p{t) [i.e. collisions) at which time one starts a new solution to (7.2) using 
the initial conditions dictated by (7.3). We also require that p > x at collision. Thus if p 
vanishes on a nontrivial interval J, we do not consider an x(t) with x(t) = p(t) for t G J a 
solution. Note that since p is sub-sine, if p > x at one collision the same inequality holds at 
the next collision, i.e. the situation in Figure 7.1 does not occur. From a more physical point 
of view, the fact that p is sub-sine implies that there are no so-called graze solutions ([BD]). 
In these solutions the particle approaches the wall, touches it for an instance matching the 
wall's position and velocity, and then pulls away without a collision. The absence of these 
solutions implies that the return time to the wall is a continuous function of x(t c ) and t c . 




Figure 7.1: Solution ruled out by sub-sine hypothesis. 



§7.2 Pictorial connection of the two systems. Fix a convex curve T in the plane with 
radius of curvature function p, and assume for simplicity that the origin is interior to T. 
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Call the orthogonal projection of T onto the x-axis Iq (see Figure 7.2(a)). The distance of 
the origin to the right endpoint of Jo is p(0) as defined in §2.1. Now imagine rotating the 
plane clockwise about the origin while keeping the direction of the x-axis fixed. If Ig is the 
projection of T onto the x-axis after we have rotated by an angle 9, then the distance from 
the origin to the right endpoint of Ig is p(0). We will think of the segments Ig as the shadow 
of r and we will rotate at unit speed and identify the angle 9 with the time t. It follows 
from (2.3) that the right endpoint of the shadow is evolving according to p + p = p. Thus 
the end of the shadow is moving like the wall in the impact oscillator. 

t <f(z) 




M ► 4 z 

x(t) 

M — 

P( 

Figure 7.2: (a)Pre-collision in the rotating plane 
model; (b) Collisions in the rotating plane model. 



To introduce the particle, pick a point zq exterior to T and monitor the evolution of its 
shadow (its orthogonal projection onto the x-axis) as we rotate the plane. The point moves 
in a circle and so if its shadow is given by x(t), then x + x = 0. 

Now for the collisions. Assume that at time t, x(t) > p(t) as in Figure 7.2(a). As we 
continue rotating, at some time t c , x(t c ) will equal p(t c ), i.e. the shadow of the point and 
that of T will collide. This will be when the supporting line to T that contains z$ is vertical. 
Since the shadows of zq and Y are taking the role of the particle and the wall in the impact 
oscillator, when the shadows collide we should act on the point zq in such a way that its 
shadow rebounds from the shadow of T with the appropriate velocity. It is clear from Figure 
7.2(b), that the correct action is to reflect zq around the point of tangency to the T to obtain 
a new point &(zo). As the plane continues to rotate, the shadow of $(^o) will evolve away 
from the shadow of T as if it had a perfectly elastic collision. The map $ just defined is 
clearly identical to the dual billiards map associated with T. 

§7.3 Connecting the generating functions. To get a more precise connection between 
dual billiards and the impact oscillator we connect the generating function of the half-cylinder 
twist map associated with dual billiards with an action integral of an impact oscillator 
equivalent to that of §7.1. 

First we derive a new formula for the generating function of dual billiards. If B(0\, 62) is 
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the area of the sector of T between the 9\ and 9 2 (see Figure 7.3), then using Stokes theorem 

B(6 1 ,e 2 ) = \ f 62 a(9) An e ds(9) 

Recalling that p = p + p", integrating by parts yields 

B(e 1} e 2 ) = 1( f e 2 p 2 - (p'f de + P (e 2 )p' \e 2 ) -p(9 1 ) p '(9 1 )). (7.4) 



B(e r e 2 ) 



Figure 7.3: The region B(9 1 ,9 2 ). 

We now move the origin to the point Q = Cg 1 n Cg 2 and let r](9) be the distance from 
Q to Cg. With this new origin the area computed as B(9\, 9 2 ) is the generating function of 
the half-cylinder twist map associated with the dual billiards map (see Remark 3.2.1). In 
a slight abuse of notation we will treat this as a function of 9 rather than x and denote it 
hr(9i,9 2 ). To use (7.4) we note that from the point of view of the point Q, the radius of 
curvature of T is —p(9), and so using the fact that 77(^1) = i](9 2 ) = 0, we get 

1 I /„/\2 J2 



h T {9i,9 2 ) = - J W-rfdB. (7.5) 

where the function rj{9) is the (unique!) solution to 

rj" + r] = -p 

with initial conditions 77(6*1) = r](9 2 ) = 0. For future reference note that from the geometric 
interpretation of the derivative of the height function (see §2.1), rj'(9i) = R(9i,9 2 ) and 
r ] '(9 2 ) = -L(9 l ,9 2 ). 

To connect this action to the impact oscillator above we define for that system a new 
variable rj(t) = x(t) — pit); thus i] measures the distance between the particle and the wall. 
The collision rule is now 

-Vbefore(t c ) = Vafteritc), (7.6) 
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and motion away from collision is governed by 



t) + r] — -p. (7.7) 

It will be useful to think of (7.6) and (7.7) as defining a new impact oscillator that is 
dynamically equivalent to the original one. In this new system, 77 is the position of a particle 
governed by (7.7) that has perfectly elastic collisions with a stationary wall. 

Now between impacts this system has Lagrangian L(x,x,t) := x 2 /2 — x 2 /2 — xp(t) and 
Hamiltonian H(q,p,t) := q 2 /2 +p 2 /2 +pp(t). If t\ < t 2 < t\ + n, we define 

hi(ti,t2) = min / L(x,x,t) dt 

where the minimum is taken over all paths x(t) with x(ti) = and x(t 2 ) = 0. Thus 

It! 



rt2 

hi(t!,t 2 ) = / L(r],f],t) 



where r\ = r)(t] t\, t 2 ) is the (unique!) solution to rj + r\ = —p with boundary values i](ti) = 
= V (t 2 ). 

Now note that hj is the Hamilton principal function (see Remark 7.3.1) with the starting 
and ending configuration variables fixed at zero, thus 

d 2 h I (t 1 ,t 2 ) = -H(r)(t 2 ),r)(t 2 ),t 2 ) = -r]{t 2 ) 2 . 

Recalling the collision rule (7.6), we see that a sequence of times (tj) will be the sequence of a 
collision times of a solution to the oscillator precisely when it is gives stationary configuration 
as in §1.2 for the variational problem defined using hj. But exactly as in (7.4) one has 

t2 r t2 2 



/•I2 ft 

/ VP = 



Hi Jt x 

and so h r {ti,t 2 ) = -hi(ti,t 2 ). 

In conclusion, given an impact oscillator described by (7.6) and (7.7), define as the 
return map to collisions, i.e. : (to, (v(to)) 2 ) i— > (ti, (v'i^l)) 2 ) where to and t\ are the times 
at consecutive collisions; we have shown that $/ is identical to the half-cylinder twist map 
map associated with dual billiards on a convex curve with radius of curvature function p. 

Remarks: 

7.3.1: If H is a Hamiltonian that satisfies the Legendre condition, then Hamilton's 
principal function is 

rh 

S(qo,qi,t ,t{) = / L(x,x,t) 

where L is the corresponding Lagrangian and x{t) is a solution to Hamilton's equations (or 
equivalently an extremal of the integral) with x(to) = go and x{t\) = q\ (see eg. chapter 9 
in [Gs]). For S to be well-defined this solution needs to be unique. The differential of the 
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principal function is dS = p dq — H dt. What was defined as hi above is the restriction 
of S to = qo — q\. In [M5] Moser uses a similar construction, but restricts S to to — 
and t\ — 1 to obtain the generating function of the time-one map of the solution flow to 
a time-periodic Hamiltonian. He shows that any twist map of the compact annulus can be 
thus obtained. This raises the question of which twist maps of the half-cylinder can arise 
from impact oscillators. 

7.3.2: One can also find solutions to the impact oscillator by considering broken paths 
consisting of a collection of functions Xi(t) and times t; so that Xi(ti) = = Xi(ti + i) and 
Xi > 0, otherwise. Solutions are then extremals for J2 L(xi,Xi,t). 

In this context introducing the generating function hj is analogous to the strategy of 
broken geodesies in the variational theory of geodesies (see page 330 of [Gl] for some history). 
For extremals of the fixed endpoint problem the introduction of hj replaces an infinite 
dimensional problem with a finite dimensional one. Crucial to the application of this method 
is the fact that that there is a unique extremal from x = to x = which begins at a time 
to and ends at ti, where to < ti < + 71 ■ 

7.3.3: Yet another way of connecting the dual billiards map to an impact oscillator can 
be obtained through a kind of regularization of the collisions with the wall. Solutions to the 
equation 

rj + rj = -sgn(x)p(t) (7.9) 

behave as if they pass through the wall instead of colliding. Upon passage through the wall 
the sign of the forcing is changed. Solutions of (7.9) for x < are the reflections through 
the origin of those for x > 0. We can use (7.9) to generate a flow on R 2 x S 1 via 

x = y 

y = -x - sgn(x)p(r) 
t = 1 

The return map of the flow to the open annulur cross section defined by x = 0, y > will 
be the second iterate dual billiards map on the T with radius of curvature p. But note that 
the map will be in (6, £) coordinates and will thus not be area-preserving. 

The correct choice of second coordinate (i.e. the one that makes the map area-preserving) 
can be seen to be dictated by (7.8) and Hamilton- Jacobi theory as in Remark 7.3.1. From 
yet another point of view we can make the 1 ^-degree of freedom system given by (7.9) into 
an autonomous 2-degree of freedom solution in the usual way by letting H*(q,p, E,t) = 
H(q,p,r) — E, where now H(q,p,t) = p 2 /2 + q 2 /2 + \q\p(t). Letting the r variable be 
on a circle, and restricting attention to H* = 0, the return map in (r, E) coordinates to 
q = 0,p > will be the second iterate of dual billiards map. Note that this map is area 
preserving because the flow preserves volume and since q = on the chosen cross section, 
the E coordinate there is just p 2 /2. 

7.3.4: In [Kl], [SV] and [GS], it is shown that dual billiards on rational polygons are 
stable, i.e. there do not exist any unbounded orbits. Recalling §2.1, this implies that a class 
of impact oscillators as defined by (7.6) and (7.7) whose forcing functions are contained in a 
certain collection of linear combinations of delta functions have no unbounded solutions. 
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7.3.5: The physical content of the fact that the return map to the wall using the (t,i] 2 ) 
coordinates is a twist map is simply that when two points are started at the same time at 
the wall, the one with greatest initial velocity will take the longest to return to the wall. 

Appendix 

The appendix is a remark on invariant circles for (inner) billiards. As is the case with 
dual billiards (§4.2), the solution to the inverse problem for the existence of invariant circles 
is a special case of Proposition 1.3. The remark rests on a elementary computation, but it 
does not seem to be widely known. Billiards is the subject of a large body of work (cf. [T3]); 
we only briefly describe the system here. Billiards as a twist map is considered in [Ml], [D] 
and [Me2]. 

If T is an oriented convex curve in the plane, the billiards map is defined using a point 
particle moving freely in the region bounded by T. The particle has perfectly elastic collisions 
with the curve, so the angle of incoming motion with the tangent to the curve is equal to 
the angle of the outgoing direction. The billiards map just keeps track of collisions with the 
curve. The map can be described using coordinates t, the arc length parameter along the 
curve, and u, the angle between an outgoing line of motion and the tangent to the curve. 
The billiards map is then (to,uo) h- > (ti,u>i), where (to,uo) represents an initial point on 
the curve and an outgoing direction, and (ti,ui) represents the point on the curve at the 
next collision and its outgoing direction after the collision (see Figure A.l). This defines a 
homeomorphism of the annulus S* 1 x [0,7r]. 




Figure A.l: The billiards system. 

An invariant circle of this annulus map corresponds to a geometric object called a caustic. 
For simplicity we restrict to the case of convex caustics. A convex caustic is a convex curve 
inside of T with the property that an outgoing path that is tangent to Ti on one pass must 
also be tangent on passes after all future (and past) collisions. The inverse problem for 
the existence of caustics is the following: given a convex curve Ti, does there exist another 
convex curve T so that Ti is a caustic for billiards on T? The solution to this problem 
involves the string body construction. Fix the curve Ti and a length of string that is longer 
than the perimeter of T\. Loop the string around T\ and put a pencil in the loop, pull the 
string tight and move the pencil around Ti keeping the string tight. The pencil will draw 
out a curve T that is characterized by the constancy of the length B + C + D as shown in 
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Figure A. 2. It is clear that T will have Ti as a caustic. The converse is also true; if T has 
a caustic Ti, then T may be obtained as a string body of Ti. Put another way, given any 
r2 inside T we can define a quantity called the string parameter L(t) = B(t) + C(t) + D(t) 
where the quantities B, C, and D may vary with the point Pit). The curve T2 is a caustic 
if and only if L{t) is constant. This result is attributed to Minasian in [S]; it also appears in 
[P] (see also [Tr]). 




Figure A. 2: The string parameter. 



With the proper choice of coordinates, the annulus homeomorphisms arising from bil- 
liards on convex curves gives a class of area preserving twist maps that have a theory similar 
to that described in §1. In particular, there is a generating function for the homeomorphism, 
and one can find orbits using variational techniques. The analog of Theorem 1.1 is true as 
well as much else. Our purpose here is just to remark on one aspect of the theory, namely, 
the connection of the area function A{t) defined in §1.4 with the string body construction. 

To describe the billiards map more carefully, begin by letting f3(i) = (Pi(t), @2{t)) be a 
parameterization of convex curve T by arc length. For simplicity, assume that T has perimeter 
2tt, and so t E S 1 . The generating function for the billiards map is the opposite of the length 
of the chord between impacts with the curve, 

hfah) = -\\p{h) - p{to)\\. 

Thus, for example, one can obtain period three orbits by finding a triangle with edges on T 
that has maximal perimeter among all such triangles. 

In accord with (1.1) we differentiate h to get the appropriate second coordinates to make 
billiards an area-preserving twist map, 

~ p{to)\\ 

where a; (to, t\) is the angle labeled ojq in Figure A.l. Thus the appropriate second coordinate 
is the opposite of the cosine of the angle between the outgoing ray and the tangent line to 
the curve. Billiards defines a twist map because h\2 < 0. The billiards on a convex curve Y 
then defines an area- preserving twist map -^r '■ S 1 x [— 1, 1] — > S 1 x [—1, 1]. 

Now assume that T has a caustic, i.e. there exists a ^r- invariant circle. By the appro- 
priate version of Birkhoff's Theorem (Theorem 1.1(b), see [Ml]), an invariant circle is the 



33 



graph of a Lipschitz function u : S 1 — > [—1,1]. The function whose graph is the lift of the 
invariant circle to M x [—1,1] will also be denoted u. The coordinates on both S 1 and R will 
be called t. 

Examining the geometric situation one sees that a caustic Ti defines a one-parameter 
family of angles each associated with a point on T. The angle at the point (3(t) is denoted 
v(t). The point of tangency to Y\ of the outgoing ray from (5{t) is denoted a{t) and the 
advance map along T is g(t) (see Figure A. 3). 




P(t) 

Figure A. 3: Functions associated with a caustic. 
Using the connection of u(t) to these geometrically defined quantities, 

u(t) = -cos(v(t)) = -P(t). Mt) _ mr 

and so 

|( Wt )- W |l) = „( t ) + 1 £|zM.v (t ). 

The area function (§1.4) associated with the invariant circle is 

rg(t) 



A(t) = / u(s) ds-h(t,g(t)) 

J t 

\a(s)-(3(s 



g(t) 

\\a'(s)\\ ds+\\a(g(t))-P(g(t))\\ 

-\wt)-p(t)\\ + \m(t))-m\\ 

where we used the fact that | | ^(g)I^(g) || is the unit vector in the direction of cr'is). 

The first term in the final expression for A(t) is the opposite of the length along T\ 
between the two tangent rays (this length is labeled E in Figure A. 2). If the perimeter of Ti 
is P, this first term is thus D — P. The second and third terms in the expression for A(t) 
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are B and C as shown in Figure A. 2, and so A(t) = L(t) — P. Thus using the appropriate 
version of Proposition 1.3, the string parameter is constant if and only if Ti is a caustic for 

r. 

Remark A.l: Note that in the billiards systems the generating function is always 
negative. The geometric interpretation of the generating function as an area in the annulus 
(Remark 1.1.1) is still valid, but now the area is below the t-axis. The area function in the 
case of billiards is a positive number giving the area below the invariant circle and between 
a vertical arc and its image. 
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